Bulk electronic structure of optimally doped Ba(Fei_ x Co x ) 2 As 2 
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We report high-resolution, bulk Compton scattering measurements unveiling the Fermi surface 
of an optimally-doped iron-arsenide superconductor, Ba(Feo.93Coo.o7)2As2. Our measurements are 
in agreement with first-principles calculations of the electronic structure, revealing both the X- 
centered electron pockets and the T-centered hole pockets. Moreover, our data are consistent with 
the strong three-dimensionality of one of these sheets that has been predicted by electronic structure 
calculations at the local-density-approximation-minimum As position. Complementary calculations 
of the noninteracting susceptibility, xo(q, cj), suggest that the broad peak that develops due to 
interband Fermi-surface nesting, and which has motivated several theories of superconductivity in 
this class of material, survives the measured three dimensionality of the Fermi surface in this family. 



PACS numbers: 71.18.+y, 74.70.-b, 78.70.Ck 

The discovery of superconductivity in LaOi-^F^FeAs, 
a member of the so-called "1111" family of iron-pnictide 
materials pQ, has triggered intense research into its na- 
ture and origin. Tempted by the proximity of the su- 
perconductivity in these compounds to an antiferromag- 
netic spin density wave (SDW) state, many theoret- 
ical models have focused on spin fluctuations as the 
key [2j [3], supported by first-principles calculations of 
the electronic structure. Such calculations predict four 
Fermi-surface (FS) sheets: two hole sheets at the center 
of the Brillouin zone (BZ), and two electron sheets cen- 
tered at its corner [4j [5]. In the undoped (and lightly 
doped) "1111" compounds the hole and electron sheets 
are well nested, leading to a broad peak in the noninter- 
acting susceptibility [2]. The spin fluctuations that are 
thought to develop from this FS instability have already 
been observed in neutron scattering measurements [6], 
and theoretical models suggest they promote a strong 
pairing interaction, provided that the order parameter 
changes phase between the hole and electron FS sheets 
(the so-called s±-model) [2j [3]. Central to such mod- 
els, however, is a detailed understanding of the topology 
of the FS; here we present a bulk measurement of the 
FS of an optimally-doped iron-arsenide superconductor, 
Ba(Feo.93Coo.o7)2As2, alongside complementary calcula- 
tions of the noninteracting susceptibility, Xo(q>^)- We 
demonstrate the data are most consistent with a three- 
dimensional (3D) FS but that the nesting that is required 
for the <s±-model survives, lending its availability to all 
iron-pnictide compounds. 

Ba(Fei_ :E Co :E )2As2 is a member of the "122" fam- 



ily of iron-pnictide superconductors, with a maximum 
T c ~ 24 K for x ~ 0.06 0|8]. Unlike the "1111" family, 
whose FS is quasi-two-dimensional (2D), the "122" com- 
pounds are predicted by electronic structure calculations 
to have substantial 3D warping of one of the hole sheets 
[5] , the degree of which is sensitive to the position of the 
As with respect to the Fe plane. Experimental studies 
have not yet adequately resolved wether the electronic 
structure is better described by the experimental As po- 
sition or that at the local density approximation (LDA) 
minimum. Nonetheless, both families exhibit nesting be- 
tween parallel pieces of the hole and electron FS sheets. 
Neutron scattering measurements reveal spin fluctuations 
at the FS nesting vector, which is coincident with the an- 
tiferromagnetic (AFM) wave vector for both "1111" and 
"122" compounds [6 . For Fe(Te,Se), where the nesting 
vector is not coincident with the AFM wave vector, spin 
fluctuations are only evident at the FS nesting vector [9] , 
indicating that the long-range magnetic order is not FS 
driven [TO] . 

Experimentally, de Haas- van Alphen (dHvA) measure- 
ments have found good agreement with LDA calcula- 
tions, observing orbits that correspond to both elec- 
tron and hole FS sheets in LaFePO [TT] and non- 
superconducting SrFe2P2 [E]- Indeed, only small rigid 
shifts of the LDA band structure were required in both 
cases to achieve good agreement with the experimental 
data. However, owing to restrictions imposed by the 
electron mean free path, these measurements have so 
far been confined to compounds away from the maxi- 
mum Tc. Angle-resolved photemission (ARPES) mea- 




FIG. 1: (Color online) (a) The LDA FS of 

Ba(Feo.93Coo.o7)2As2 (top) and its relation to the occu- 
pation density projected down the c* axis (bottom). Note 
that the occupation density has been convoluted with the 
experimental resolution function. (b) The overlaid FS 
contours (black and white lines, top) are identified via the 
maxima in (the absolute value of) the first derivative of the 
occupation density shown at the bottom along a projected, 
high symmetry path. 

surements, on the other hand, yield a more contradictory 
picture. Whereas some measurements on the undoped or 
lightly-doped BaFe2As2 system demonstrate reasonable 
broad, albeit qualitative, agreement with LDA predic- 
tions [T3HT5] , others suggest either the existence of hole- 
like "blades" in place of the quasi-2D electron sheets ex- 
pected at the corner of the BZ [16] that are interpreted 
as indicative of additional order, or do not simultane- 
ously observe both hole and electron FS sheets. Indeed, 
scanning-tunneling microscopy has been applied to the 
surface layers of Ba(Feo.93Coo.o7)2As2, finding a complex 
array of sample-dependent (specifically, cleavage temper- 
ature) ordered and disordered structures [17], which are 
not representative of the bulk electronic structure. 

More recently, several ARPES studies [I3j [I8j [19] , sup- 
ported by dHvA [12] measurements, have reported on the 
3D nature of the "122" family FS, and particularly the 
warping of one of the hole FS sheets that has been pre- 
dicted by electronic structure calculations [5 . The FS 
that has motivated spin-fluctuation models of supercon- 
ductivity has so far been quite 2D, employing either the 
quasi-2D band structure of LaFeAsOi-^F^ [2 or a 2D 
minimal model of it [3 . Since the warping is expected 
to extend over some quarter of the c* axis, the question 
of whether the nesting instabilities that are expected to 
drive the spin fluctuations are robust against such three- 
dimensionality must be raised. 

Compton scattering provides a valuable complement 
to dHvA and ARPES, unambiguously probing the occu- 
pied states of the bulk electronic structure [20] . Its insen- 
sitivity to defects or disorder make it an ideal technique 
for tackling the FS of disordered systems (e.g. [21 ) and, 
since it is sensitive to the bulk 3D momentum density, it 
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FIG. 2: (Color online) The occupation density of 
Ba(Feo.93Coo.o7)2As2 from (a) experimental Compton scat- 
tering measurements at 295 K and (b) the electronic structure 
calculations of Fig.JIJi from the original calculation (top), and 
after rigidly shifting the bands (see text, bottom). The first 
derivative of the experimental data is shown in (c) in the same 
way as in Fig. [Ij). In (d) the FS corresponding to the shifted 
bands is shown for comparison with Fig. [TJi. 



is able to provide information about the FS throughout 
the BZ. In this study the occupation density, featuring 
the projected Fermi breaks, serves as the key quantity 
for assessing the FS of Ba(Feo.93Coo.o7)2As2. 

A Compton profile represents a double integral [one- 
dimensional (ID) projection] of the underlying 3D elec- 
tron momentum density (EMD) in which the FS mani- 
fests itself via discontinuous breaks. Six Compton pro- 
files equally spaced between [100] and [110] were mea- 
sured at room temperature (T = 295 K) on the high- 
resolution Compton spectrometer of beamline BL08W 
at the SPring-8 synchrotron [22], with a resolution full 
width at half maximum at the Compton peak of 0.105 
a.u.(l/8th of the BZ). For each Compton profile, ~ 
480 000 counts were accumulated in the peak data chan- 
nel, and each Compton profile was corrected for possi- 
ble multiple-scattering contributions [23]. Details of the 
sample growth and characterization may be found in Ref. 

For the purpose of comparison with the occupation 
density predicted by band theory and to aid interpre- 
tation of our data, we performed ab initio calculations 
of the electronic structure of Ba(Feo.93Coo.o7)2As2 us- 
ing the full-potential linearized augmented plane wave 
method [24] within the LDA, where the effect of dop- 
ing was included via the virtual crystal approximation 
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(VCA) [26]. The As positions are those corresponding 
to the LDA minimum [5] . Fig. |TJl shows the calculated 
FS, composed of two quasi-2D electron sheets and two 
hole sheets, in agreement with previous calculations for 
similar Co concentrations [5j [27] . Note the considerable 
three-dimensionality ("warping") of the outer hole sheet. 

Using tomographic techniques [28] , it is possible to re- 
construct the 2D EMD from a specially chosen set of ID 
Compton profiles. Here, our six Compton profiles were 
reconstructed to obtain the 2D EMD projected along the 
c* axis, which was then folded back into the first BZ [29 . 
This culminates in a projected occupation density, shown 
for the data in Fig.[2^i, representing a projection of the oc- 
cupied electron states in the BZ, where light (dark) colors 
represent high (low) occupation densities. As can be seen 
from the data, both the electron and hole FS sheets that 
are predicted by electronic structure calculations (Fig. 
^p) are clearly identifiable, with the electron sheets in 
the corner of the projected BZ (XP) and the hole sheets 
in the center, YZ. The slight flattening in the measured 
occupation density midway between YZ and XP arises as 
a consequence of describing its rapidly varying anisotropy 
with a finite number (six) of measured profiles and does 
not interfere with the location of the FS breaks. It is 
worth emphasizing that although this quantity is 2D, in- 
formation on the third, projected dimension is still stored 
in the integral. For an ideally 2D FS, the Fermi breaks 
would be realised as a step function between high and low 
occupation densities broadened only by the experimen- 
tal resolution; however, the three-dimensionality of the 
FS results in an additional smearing of these projected 
breaks. 

In order to refine this picture, a method for detecting 
the FS breaks via the first derivative of the occupation 
density, as illustrated in Fig.[TJ) for the theoretical result, 
was employed. The subtraction of an isotropic compo- 
nent, corresponding to the tightly-bound electron states, 
from the data before applying the folding procedure [29] 
yields further enhancement in locating the Fermi breaks 
and hence the maxima in the first derivative of this quan- 
tity were used to locate the breaks (see Fig. [2^), and the 
occupation density was contoured at the corresponding 
level [30 . Fig. |2^i shows overlaid on the data the exper- 
imental Fermi-surface breaks mapped out utilizing this 
method. Owing to the experimental resolution and the 
projected nature of this result, two features should be 
noted: first, the two electron (hole) sheets are not re- 
solved individually and second the contours here repre- 
sent an average over the shape and size of the FS sheets 
projected down the c* axis. The slightly larger size of the 
electron compared with the hole pockets is expected as a 
consequence of electron doping. For these (experimental) 
contours assuming a quasi-2D FS topology, this leads to 
an occupied volume of the BZ of 4.38 ±0.16 electrons. In 
comparison, for Ba(Feo.93Coo.o7)2As2, the FS should en- 
close 4.14 electrons, and we now focus on an explanation 



for the discrepancy between these two figures. 

In Fig. [2)3, the same procedure that was applied to the 
data has been employed for the theoretical FS, which 
explicitly includes the three-dimensionality of the outer 
hole sheet, and good agreement is found between exper- 
iment and theory. However, our method for locating the 
Fermi breaks relies on a peak in the first derivative of 
the occupation density, and as demonstrated in Fig. [TJd, 
no such peak is discernible in the theoretical occupa- 
tion density that corresponds to the wider part of the 
"warped" hole FS. In projection, and with our experi- 
mental resolution, this feature is not resolvable and is 
therefore expected to be absent in the experimental con- 
tours as well. On the other hand, since some of the hole 
FS is not accounted for in the contours, the "warping" 
still leaves its mark in the occupation density in the form 
of an imbalance in the detected electron count. Indeed, 
we find that 4.34 electrons are enclosed by our theoret- 
ical contours, where the excess electrons stem from an 
underestimation of the volume of the hole FS, i.e., as a 
direct consequence of the "warping" of one of those cylin- 
ders. Since the theoretical 3D FS used in this analysis 
unambiguously encloses 4.14 electrons, this supports well 
the notion that the extra electrons in our experimental 
contours provide strong evidence of the 3D warping of 
the hole FS. It should be noted that the theoretical FS is 
calculated with the As atoms at the LDA minimum [5]; 
crucially, our calculations at the experimental As posi- 
tion, in which the "warping" is much weaker, were not 
able to explain such a mismatch in the occupied fraction 
of the BZ. 

Turning our attention to the details of the FS con- 
tours, and their comparison with theory, our data sug- 
gest slightly larger FS sheets than those expected from 
the calculation. It is instructive to investigate whether 
such a discrepancy can be accounted for by a rigid shift 
of the theoretical band energies, in such a way that the 
electron sheets increase in size while the overall electron 
count in the system is maintained by a shift on the hole 
sheets. As can be seen in Fig. [2}d, shifts of ~ — 120 me V 
and ~ +80 meV for the electron and hole bands, respec- 
tively, lead to a good overall agreement between data and 
theory in which the shifted theoretical contours now en- 
close 4.40 electrons. As illustrated in Fig. [2]i employing 
such shifts preserves the FS topology, only affecting the 
size of the hole and electron pockets. 

Overall, the projected FS topology obtained in our 
study is in good agreement with other measurements on 
similar compounds (e.g. p~2j [13] US]), as well as elec- 
tronic structure calculations within the LDA. In Ref. 
P2], a rigid shift of the theoretical electronic structure 
was required to well reproduce their quantum oscillations 
measurements on SrFe2P2- Those shifts were of a simi- 
lar magnitude, and were in the opposite direction; how- 
ever, the necessary additional approximations employed 
in this calculation, specifically the VCA which is an aver- 
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FIG. 3: (Color online) The real and imaginary parts of the 
static susceptibility, x(q 5 0)? of Ba(Feo.93Coo.o7)2As2. The 
wave vector q is expressed in units of 2n/a. 



age potential approximation may not provide a complete 
picture of the disorder in contrast to a more sophisti- 
cated approach such as the coherent-potential approxi- 
mation [31]. Additionally, our data show no evidence of 
the propeller-like structure near the X-point that has re- 
cently been reported in some ARPES measurements [16] . 

The nesting between these hole and electron FS sheets 
has been implicated in several models of superconductiv- 
ity, motivated by a quasi-2D FS topology that provides 
ample nesting throughout the c* axis of the BZ [2j[3], as 
evidenced by the broad peak in the susceptibility near 
q = (7r,7r) (2j. However, for the "112" compounds pre- 
sented here, our data, in agreement with some ARPES 
and dHvA studies p~2j [13] [18] [19] , suggest that an appre- 
ciable portion of the FS along the c* axis does not nest at 
the same wave vector. In order to address what impact 
this may have on the susceptibility of the system, we have 
calculated the noninteracting susceptibility, Xo(q>^)> f° r 
wave vector q and frequency uj [2 j . A previous calculation 
of Xo(q> kj) of the "122" materials [32] has focussed on the 
stoichiometric and hole-doped compounds, and with the 
As atom located at the experimental position, for which 
the three-dimensionality of the Fermi surface is substan- 
tially less pronounced. Both real and imaginary parts 
of xo(q^) were calculated for the (raw, unshifted) elec- 
tronic structure of Ba(Feo.93Coo.o7)2As2, and are shown 



in Fig. [3] Despite the three-dimensionality of the FS, a 
substantial (although broad) peak in the imaginary part 
of Xo(q) 5 whose origin lies in FS nesting, survives into the 
real part, with a maximum at q max ~ (0.56, 0.56)7r/a. 
This broadly peaked structure bears close resemblance 
to that calculated for the "1111" compounds in Ref. [2], 
although the strength of the peak is slightly weaker here 
due to the "warping" in k z . The wave vector of the 
maximum itself, q max agrees favorably with the nest- 
ing vector obtained from our experimental contours of 
Qexp ~ (0.62, 0.62)7r/a, given that our contours repre- 
sent an average of the FS sheets [applying the same pro- 
cedure to the convoluted theoretical occupation density 
yields q t h e ~ (0.61, 0.61)7r/a]. That our experimental re- 
sults are in good agreement with the LDA calculation, 
particularly with the topology of the FS, provides exper- 
imental evidence that the s±-model may still provide a 
route to superconductivity in the "122" compounds, de- 
spite the enhanced three-dimensionality of the electronic 
structure. 

In conclusion, the FS topology of Ba(Feo. 93000.07)2 As2 
has been investigated using Compton scattering, a true 
bulk probe of the occupied electron states. Good agree- 
ment is found between our experimental projected FS 
and the predictions of the LDA, which can be max- 
imised by employing only small rigid shifts of the the- 
oretical electronic structure. In particular, our data 
cannot be explained without including the strong three- 
dimensionality of the hole FS expected for the "122" fam- 
ily at the LDA-minimum As position (rather than the ex- 
perimental one). However, our calculations of the nonin- 
teracting susceptibility of this compound suggest that the 
peak that has been implicated in some theories of super- 
conductivity survives this three-dimensionality and may 
still be available as a mechanism in these compounds. 
Such an observation provides experimental evidence that 
the s±-model (and others that rely on FS nesting) may 
be general models for superconductivity in the iron pnic- 
tides. 
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